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The use of quasi-elastic electron nucleus scattering is shown to provide significant 
constraints on models of the proton electromagnetic form factor of off-shell nucleons. 
Such models can be constructed to be consistent with constraints from current con- 
servation and low-energy theorems, while also providing a contribution to the Lamb 
shift that might potentially resolve the proton radius puzzle in muonic hydrogen. 
However, observations of quasi-elastic scattering limit the overall strength of the 
off-shell form factors to values that correspond to small contributions to the Lamb 
shift. 
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I. INTRODUCTION 

The structure of nucleons bound in the nucleus is different than that of free nucleons. A 
prominent example is the EMC effect in which the influence of the medium is known to modify 
the quark distribution functions of nucleons, see the reviews [T]. An even older example of 
medium effects involves the neutron which lives forever in stable nuclei. This is because the 
effects of the pe~v component of the free neutron wave function are suppressed by the effects of 
binding energy in nuclei. The existence of medium effects on the structure of nucleons cannot 
be denied, but elucidating all of the possible effects and the relations between them is a task 
for ongoing research. We shall focus here on medium modifications of proton electromagnetic 
form factors [2], which potentially affect quasi-elastic scattering [3], [I], and may contribute 
to solving the proton radius puzzle [5]. 

A prominent example of medium modifications is the work of Ref. [2] which involves mea- 
suring the double ratio of proton-recoil polarization-transfer coefficients of the quasielastic 
4 He(e,e'p) 3 H reaction with respect to the elastic 1 H(e,e'p) reaction which is sensitive to pos- 
sible medium modifications of the proton form factor in 4 He. Measurements of this double 
ratio at four-momentum transfers squared between 0.4 GeV 2 and 2.6 GeV 2 performed at both 
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Mainz and Jefferson Lab find a reduction of about 10% in the double ratio, which corresponds 
to a similar reduction in the ratio of electric to magnetic form factors Ge/Gm- Models which 
treat the nucleon as a bound state of three quarks which move under the influence of quarks 
in other nucleons, consistent with the EMC effect and much nuclear phenomenology, predict 
such a reduction |6J. Alternative explanations involving final state interactions are possible 
(see the discussion in Ref. [2] and references therein), but seem to be incompatible with the 
totality of relevant data. 

It is noteworthy that measurements of quasi-elastic scattering are related to one of the most 
perplexing physics issues of recent times-the proton radius puzzle. This puzzle originates in 
the extremely precise extraction of the proton radius [7] from the measured energy difference 
between the 2Pf^ 2 and 2S[^ 1 states of muonic hydrogen (H). This Lamb shift depends on 
the finite size of the proton's electromagnetic field. The extreme precision of the measurement 
leads to an extracted value of the proton radius smaller than the CODATA |8J value (extracted 
mainly from electronic H) by about 4% or 5.0 standard deviations. This implies [7J that 
either the Rydberg constant has to be shifted by 4.9 standard deviations or that present 
QED calculations for hydrogen are insufficient. Since the Rydberg constant is extremely 
well measured, and the QED calculations seem to be very extensive and highly accurate, the 
muonic H finding has presented a significant puzzle to the entire physics community. 

We need a brief discussion of the relevant phenomenology. Pohl et al. show, perturbatively 
that the energy difference between the 2P 3 ^ 2 and 2S[^ states, AE is given by 

AE = 209.9779(49) - 5.2262r 2 + 0.0347rJ meV, (1) 

where r p {r 2 is related to the slope of Ge{Q 2 ) at Q 2 = 0) is given in units of fm. Using this 

equation, and the experimentally measured value, AE = 206.2949 meV, one can see that the 
difference between the Pohl and CODATA values of the proton radius would be removed by 
an increase of the first term on the rhs of Eq. (1) by just 0.31 meV=3.1 x 10~ 10 MeV, but an 
effect of even half that much would be large enough to dissipate the puzzle. 

This proton radius puzzle has been attacked from many different directions [9] [22], [5] The 
present communication is intended to investigate the hypothesis [5j that the off-mass-shell 
dependence of the proton electromagnetic form factor that occurs in the lepton-proton two- 
photon exchange interaction can account for the 0.31 meV. This idea is attractive because the 
computed effect is proportional to the lepton mass to the fourth power, and so is capable of 
being relevant for muonic atoms, but irrelevant for electronic atoms. To make a calculation 
one needs to postulate a specific dependence of the electromagnetic form factor as a func- 
tion of the proton's virtuality (difference between the square of the proton's four-momentum 
vector and the square of the proton mass). Many different functional forms are possible 
and Ref. |5] chose one that accounted for the difference between the muonic and electronic 
hydrogen measurements. 

Ref. [21] uses a dispersion analysis of the two-photon exchange term [25], [26] to provide 
limits on the size of the allowed off-shell effect of the specific chosen form of Ref. [3] . This is 
done by expressing the virtual Compton scattering amplitude implied in Ref. [5] in terms of the 
invariant and relating those amplitudes to electric and magnetic polarizabilities. The ones 
used in Ref. [5] are shown in Ref. [21] to be far larger than the experimentally measured ones. 
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The accuracy of such dispersion relation approaches may be less than previously thought [27] . 
Nevertheless, we construct a new model of the off-shell form factor that is consistent with all 
of the conditions mentioned in Ref. [21]. These conditions are derived using the constraints 
of second order in chiral perturbation theory. Very recently Birse & McGovern [2H] evaluated 
the constraints to fourth-order in chiral perturbation theory. We show how to develop off-shell 
models that are consistent with any order of chiral perturbation theory. Moreover, these kinds 
of models are testable in a variety of arenas, and in particular it is of interest to examine the 
consequences of the proposed off-shell model of [5] for electron-nucleus scattering. 

The idea we consider is that a bound nucleon can be taken evanescently off its mass- 
shell by virtue of its interactions with other nucleons, and that the consequences of using 
the model [5] can therefore be tested. We make an explicit calculation of how the ratio of 
proton electromagnetic form factors Ge/Gm is changed in the nuclear medium according 
to the model of Ref. [5], and confront the ensuing predictions with the data of Ref. [2], 
Sect. II. The model of [5] is shown to yield medium modifications of the ratio of Ge/Gm far 
in excess of the observed effects. A new model is constructed in Sect. Ill that leaves the ratio 
Ge/Gm unmodified in the medium. This model also is constructed to be consistent with the 
restrictions of any finite order in chiral perturbation theory. The off-shell modification depends 
on an overall strength parameter denoted as A, which is limited by quasi-elastic scattering. 
The corresponding change in the Lamb shift is computed in Sect. IV. We find that the use of 
the largest values of A allowed by quasi-elastic scattering lead to changes in the Lamb shift 
that are far too small to account for the proton radius puzzle. 



II. OFF-SHELL PROTON FORM FACTOR IN QUASI-ELASTIC ELECTRON 

SCATTERING 

The version of the Dirac form of the electromagnetic vertex operator for an interaction 
between one on-mass-shell and one off mass-shell nucleon used in Ref. [5] can be expressed as 



(p + p'Y f s - M 



1 



= 7" Fife' 2 ) + ST" (3) 
2M M K J 



-\ q 
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for a photon interacting with an proton that is initially off its mass shell, where M is the 
nucleon mass, p' = p + q, either p^ or p /fl are off the mass shell and p oS is the four momentum 
of the off shell nucleon ((p oS ) 2 ^ M 2 ). Note that this was one of the three possible forms (O-i) 
of operators proposed in Ref. [5]. Other forms are possible. The values 



(6) 
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were used in [5] to give a contribution to the Lamb shift large enough to allow the CODATA 
value of the proton radius to be consistent with the Pohl experiment. We note that the use 
of Eq. ^ in [5] was consistent with current conservation. Replacing (p + by (p + p')^ — 
(p + p') ■ qq^/q 2 (where q is the virtual photon momentum gives no change to the computed 
shift in the atomic binding energy. 

The use of the vertex function of Eq. ^ in computing virtual-photon-proton Compton 
scattering leads to new contributions at low values Q 2 = — q 2 . It is worthwhile to compare 
these effects with those of standard formulations in which two invariant amplitudes T 12 appear. 
Given the model of Eq. ^ there is a new contribution to T2 (but not Ti), which has been 
found to be Ref. [23] 

for small values of Q 2 . In standard formulations the coefficient of the Q 2 term of T2 is given 
in terms of the electric and magnetic (3m polarizability of the proton as Q 2 /e 2 (aE + (3m)- 
Equating the coefficients of Q 2 gives A/6 2 = — 0.018(2)/(79 MeV) 2 . Thus the constraints 
imposed by the known electromagnetic polarizabilities yield a value of A/6 2 that is about 55 
smaller than and of the opposite sign to the value given in Eq. (JgJ) [24J. The other models 
mentioned in Ref. [5] were not used to compute the Lamb shift, but would correspond to 
different values of A/6 2 which are of the same order of magnitude as that of Eq. and 
those models would therefore fare equally poorly. Thus the model of Ref. [5] is not consistent 
with known features of the virtual-photon-proton Compton scattering amplitude. However,it 
is worthwhile to examine the consequences of such a model for other processes to illustrate 
the connections between different areas of physics. Moreover, the restrictions of Ref. [23] can 
be removed simply by postulating that the off-shell effects of Eq. §2§ be proportional to q 4 . 

We therefore consider quasi-elastic electron-nuclear scattering. The basic idea is that a 
proton, bound in the nucleus is slightly off its mass shell. The average binding energy of a 
nucleon is less than 1 percent of its mass. But even this very small binding corresponds to an 
off-shell effect that contradicts experiment if the model of Eq. ^ is used. Consider a proton 
bound via a Dirac mean- field Hamiltonian. The bound-state wave function obeys the 
equation 

W) (8) 



i>- m 



as in Fig. 1, in which the presence of the residual nucleus is represented by the interaction V. 
The change in the scattering amplitude dM.^ caused by the off-shell term of Eq. ^ is given 
by 

5M» = W\5T*\i/>) = F ^'\ { ^^^), (9) 

in which the final state single-particle wave function is represented by \ip')- 

We explain the relationship between the amplitude of Eq. ^ (inherent in Fig. 1) and the 
analysis of Ref. [5]. There are two interaction vertices in Fig. 1, one involving the photon and 
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FIG. 1: (Color online) A photon (wiggly line) interacts with a bound nucleon that is off its mass 
shell because of the interaction V. 



one involving the strong interaction field. The intermediate virtual proton propagator appears 
between these two vertices. This propagator is cancelled by the inverse propagator in the term 
5T, so that effectively one sees a contact interaction between the virtual photon, the struck 
proton and the residual nucleus. The interaction between the photon and the virtual proton 
converts the off-shell proton into its final state, \if)')- Thus the very same off-shell interaction 
of Ref. [5] enters here. Here one sees the combination of ST^ and V; in the two-photon 
exchange term we use a combination of 5V 1 and T u . 

The evaluation of the consequences of using Eq. (|2j) and Eq. ^ proceeds by first examining 
the single particle wave functions We treat the final state wave function as a plane 

wave represented by a Dirac spinor. This is reasonable for the present purpose because the 
effects of final state interactions are included and removed in the experimental analyses [2]. 
We aim to consider quasi-elastic electron scattering in the kinematic regime in which the 
impulse approximation is valid, and thus use a relativistic Fermi gas model to approximate 
the initial nuclear wave function. This often used approximation [29] is accurate enough for 
the schematic estimate that is the present aim. This is because any correction terms are of 
order V/M of the terms we evaluate. Under the stated assumptions we find 

( v _l V 'Y y 

5M» « F^PT^i^M^- (10) 

The largest effects of the interaction V occur at the center of the nucleus, where the density 
and the mean field can be regarded as constants. Thus we take V to be number (representing 
an average nuclear interaction) not an operator, define e = X and use — q 2 = Q 2 > to obtain 

6M» « fMp) ( {l i^f(Q 2 )) «(P)> (u) 



where 



f AQ 2 

f{Q 2 ) = ; , ^ /A9 - (12) 



1 + Q 2 /A 2 ' 
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Under the stated approximations, the present calculation is consistent with current conserva- 
tion. Replacing (p + p')^ by (p + p')^ — (p + p') • qq^/q 2 gives no change to the matrix element 
of T M appearing in Eq. ( 10 ) because the operator is evaluated between on shell spinors so that 
(p + p') 



Q 



0. 



We gain insight by using the Gordon identity to make the replacement: 
t^, so that 



{p+p'Y 

2M 



2M 



5M» = F 1 u(p') 



7 "(i + /(g 2 ))-* 



o^ v q x 
2M 



f{Q 2 



u{p), 



which shows that the nuclear medium modifies both F\ and F 2 : 

^(g 2 ) = F 1 (g 2 )/(g 2 ), <5F 2 (g 2 ) = -F 1 (g 2 )/(g 2 ) 



(13) 



(14) 



so that the change in F\ is the negative of the change in F 2 . 

We aim to see whether such modifications are consistent with present observations. Strauch 
et al. [2] measured the ratio of polarization transfer in the 4 He nucleus to that of a nucleon 
for 0.4 < g 2 < 2.6 GeV 2 . They observed a decrease of about 10%. If final state interactions 
are properly accounted for, this is a measurement of how the ratio Ge/Gm is changed in the 
medium. We therefore study the variation of that ratio. Recall the definitions 



G, 



Ft 



AM 2 



F2; G M — Fx + F 2 . 



(15) 



The medium modified form factors Ge,m are given by adding the changes in F 12 indicated by 
Eq. Q and Eq. 0. Note that G M = G M . 



The medium modified ratio is given by 



G E G E + FJ(1 + 



G 



AM' 2 ■ 



G, 



G 



G 



M 



1 + A/a + JTL- 

G E K 4M 2 ' 



(16) 



We now evaluate the function /. Our aim is to see if the smallest possible values of / are 
consistent with observations. Therefore we take e to be the ratio of the average nuclear binding 
divided by the nucleon mass (7 MeV for 4 He), so e w -0.007. Using Eq. @ we find 



f(Q 2 



-11 



A 2 



g 2 /A 2 



(17) 



which ranges between -0.6 and -1.3 as Q 2 varies between 0.4 and 2.6 GeV 2 . This is between 
6 and 25 times the effect observed by [2], if one asserts that the entire 10% reduction of 
the double ratio of polarization observables is a true medium modification. Otherwise, the 
discrepancy would be larger. 

One could argue that the model used to evaluate the nuclear effect, taking V/M to be 
a constant, is too simple to be used. The interaction V represents the nuclear mean field 
and has a spatial extent corresponding to the size of the entire nucleus. Treating this as a 
constant means that we are computing form factor modifications near the center of the nucleus. 
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This is appropriate because the experimental analyses of [2] to which we compare includes 
corrections so as to approximate the situation near the center of the nucleus. The most 
evident improvement would allow V to have an attractive scalar term and a repulsive vector 
term. Using this would lead to a larger computed effect because the cancellation between 
these terms that lead to the small average binding energy of 7 MeV per nucleon would be 
somewhat disrupted. Using V/M = —0.007 minimizes the size of the effect of using Eq. (J5]) 
in the nuclear medium. Even with this minimization, the predicted modification of the ratio 
of electric to magnetic form factors is too large. 

Note that the modified ratios that we compute do not show up in full strength in the 
(e,e'p) experiment. This is because the reaction may occur at the edge of the nucleus. But 
such effects are far too small to account for the order of magnitude problems we encounter. 

The model embodied in Eq. ^ can be regarded as ruled out by the data of Ref. [2]. The 
next section is concerned with deriving a new model. 



III. NEW MODELS THAT DO NOT MODIFY RATIOS OF FORM FACTORS 



An alternate approach is to consider the Strauch data to be a constraint. Then we have 



G, 



G 



M 



G 



M 



(18) 



where the approximation means within about 10%. We express this in terms of iq^, SFi^, 
with Fi = Fi + SFi where 5Fi being the change in Fi induced by the medium. Using the 



definitions, Eq. (15), allows us to re-express the constraint Eq. (18) as 



SJ2 

F 2 ' 



The medium modification of the ratio F 2 /Fi is experimentally accessible 
Eq. (pi leads to 



F 2 + SF 2 
F x + SFx 



F 2 
Fx 



(19) 

. The use of 
(20) 



The results Eq. (19), Eq. (20) show why our medium modification is so large. Eq. (14) 
shows that 8F\ = —5F 2 . 

The next step is to see if one can construct a model of off-shell form factors that satisfies 



O d = 



the constraints of Eq. (18)-Eq. (20). This can be done if we include an effect that changes 



F 2 so that Eq. (19) is satisfied. We can do this by adding a new off-shell term of the form 
'■" s — M) -• -. In particular, we postulate a new version of the off-shell vertex 



2M 



intended to replace the ruled-out term 5T fl of Eq. Defining this operator as we try 



0» = AF(Q 2 )[F 1 (g 2 )( 7 ^ t 



it 



) + i 



a 



2M 



''-F 2 (Q 2 )1 ^ 



'off 



M) 



M 



(21) 
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The aim is simply to find an off-shell modification that satisfies all of the constraints. Current 
conservation is explicitly satisfied by both terms. When one calculates the diagram of Fig. 1, 



the term gf q^/q 2 does not contribute because it acts between u(p + q) and u(p), see Eq. (10). 



With Eq. ([21]) we have 

8F 1 = XFF l7 5F 2 = XFF 2 , (22) 



so that Eq. (19) is satisfied. 

The use of this model in the diagram of Fig. [T] leads to an extremely simple evaluation of 
the modified quasi-elastic cross section. The effect of the medium modification is to simply 
multiply the computed quasi-elastic scattering cross section by a factor of (1 + eXF(Q 2 )) 2 . We 
thus are able to obtain a constraint on the product eXF(Q 2 ) without specifying any of the 
individual factors. The form of F(Q 2 ) is needed to compute the contribution to the Lamb 
shift and is discussed below. 

If we assume that a 10% change in the cross section (which is the typical uncertainty in 
the computation of a cross section) can be tolerated, we find that \eX\F(Q 2 ) < 0.05, or 

|A|F(Q 2 ) < 7, (23) 

for Q 2 < 10 GeV 2 . Quasi-elastic experiments have not been performed for larger values of Q 2 . 
We take e = —0.007 to obtain the above constraint. We note that this number is the smallest 
conceivable magnitude that one could extract from nuclear physics. A more detailed analysis 
would lead to a number that is larger in magnitude, and an even stronger constraint on the 



value of A. However, the limit Eq. (23) leads to a very small contribution to the Lamb shift. 



IV. LAMB SHIFT CALCULATION 



by 



The invariant lepton-proton scattering amplitude arising from two photon exchange is given 



M 



(27T) 



L^{k)T^{k,p) 



(k 2 + te) 2 



(24) 



where p is the proton momentum, and is evaluated in the common rest frame. The spin- 
averaged lepton tensor L^ v is given by 



L ^ = ^W +m \k 2 -2l.k + ie) 



(25) 



where m is the lepton mass. The term T^ v is the virtual photon nucleon forward scattering 
amplitude. We use the definition 



T MM = _ h ^-)T x {k\ k 2 ) + ^ - ^*")(p" - P ^)T 2 {k^ k 2 ) 



with k° = k ■ p/M. Then the use of Eq. (25) in Eq. (24) leads to the result: 
2m e 4 



M 



(2vr) 4 J (k 4 - (21 ■ k) 2 )(k 2 + ie) 2 



[-(2k 2 + k 2 )T 1 + (k 2 -k 2 )T 2 } 



(26) 



(27) 
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We shall be concerned with the change in T^ v caused by off-shell form factors, and denote 
the corresponding contribution to the Lamb shift, AE Q s. We use the standard procedure in 
which the zero-energy, constant amplitude M. is treated as a delta function at the origin in 
coordinate space as so that 

. d (m r a) 3 , , 

AE oS = -iM oS ^-^, (28) 

where the factor appearing to the right of Ai s is the square of the 2S muonic hydrogen 
wave function at the origin. The change in Compton scattering by our postulated off-shell 
effects is obtained by computing Compton scattering in the Born approximation. Define the 
conventional electromagnetic vertex operator for the absorption of a photon of momentum k 
as T^(k). Then 

= TZ + T G 7 = (29) 

= (^(-k) + i— (r (k) + O v {k)) + \pL -+ v, v -+ /x, k^ -k\. (30) 

= TZ + {T»{-k)O u {k) + 0"(-k)T v (k) + T u {k)0»{-k) + O v (k)T ll (-k)) 

+ [0"{-k){^+ ft + M)O v {k) + O v {k)(i>- 1/6 + M)O^(-k)} (31) 

We need the spin average, obtained by multiplying the above by ( p 1 + M) / (4M) and taking 
the trace. In the following text all quantities X^ u are spin averaged. 
We proceed by breaking up T^ u as a sum of three terms: 

= TZ + 5T^ + SZ^, (32) 

respectively of order A , A 1 and A 2 . Then we find 

rp _ F 1 2 (2p-k) 2 + 2F l F 2 k* + F 2 2 k 4 
lon ~" " M{k±-{2p-kY) " ' 1 ' 

-AF 2 M 2 k 2 + F 2 k* 
2on M{{2p-kf-k 4 ) ' 1 } 

which are standard results, and 



5T^ = Tr 



We find 



4M 



(r"(-fc)<T(fc) + 0»{-k)T v {k) + T v (k)0»{-k) + O l/ {k)Y ti (-k)) 



(35) 



ST ~AF(-fc 2 ) F 2 k 2 + AF 2 M 2 

5T i = U 77* > 5T 2 = °> ( 36 ) 



M M 2 
The second-order terms are obtained to be 



5Z »v = Tr 



(0> i (-k)(j+ ^ + M)O u {k) + O u {k){p'- }6 + M)0»{-k)) 



AM 

S7 x2p2 F 2 {F 2 {k 2 -k 2 )-2F l k 2 ) x2p2 ^F 2 M 2 -F 2 k 2 ) 
5Zi = X F — , 5Z 2 = A F — 3 . (38) 



(37) 
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The low-energy theorem and constraints of chiral perturbation theory constrain Ti(u,Q 2 ) 
for small values of v and Q 2 . Those constraints, as applied in Ref. [21] and earlier works, are 
not modified if we choose F(—k 2 ) ~ fc 4 for small values of k 2 . Thus we use 



F(-k 2 



(k 2 /A 



2\2 



;i + (-p)/a 2 )) 2 . 



(39) 



Birse & McGovern [2H] have provided constraints to fourth-order in chiral perturbation theory. 
In general, one can satisfy the constraints to n'th order by using a more general version of 



F(-k 2 ), F n (-k 2 y. 



FJ-k 2 



-k 2 /A 



2\n 



(1 + (~k 2 )/A 2 )) n ' 

Now evaluate the integral by Wick rotation 

k iK , k^ K, k 2 -)■ -K 2 - K 2 = -K 2 , K = Kcosip, \K\ = K simp 



Integrate on ip from to 7r, 



(Pk > 4vri / dK K 



<h\) sin 2 ip 



(40) 



(41) 



(42) 



Use e 2 = 47ra so 



= M^ smn I dK K 
(2tt)4 



dip 



sin 2 if) 



K 4 + Am 2 K 2 cos 2 ip 



[Ti (2 cos 2 ip + 1) - T 2 sin 2 ip)} (43) 



Now use Eq. (36) and Eq. (38) in the above to get the off-shell correction. We need 
STl + SZl = -^ F (K 2 ) -F 2 K 2 + AF 2 M 2 + ^^F, (F 2 (K 2 sin 2 ^) + 2F^ 2 ) 



M 



M 2 



M 3 



(AF 2 M 2 + F 2 K 2 ^ 
W 



5T 2 + 5Z 2 = 5Z 2 = X 2 F 2 



(44) 



8a 2 C r sin ib 

5M oS = i—mj dK K j o diP RA - 4m J 2 cog2 - m + 5 Zl )(2 cos 2 ^ + 1) + 5Z 2 sin 2 V>)] 

(45) 



The above result, along with Eq. (28), determines the value of the computed energy shift 



arising from the off-shell effect. The principal parameter is A, constrained to be less than about 



7 (Eq. (23)) from quasi-elastic scattering data. The proton electromagnetic form factors Fi j2 



are taken as dipole forms with A = 0.841 GeV, and F 2 (0) = 1.79. We start by using Eq. (39) 
and display numerical results for values of A between and 200 are shown in Fig. [2} With 
A = 7, we obtain a shift of -0.001 meV, which is about 100 times too small to significantly 
affect the Lamb shift calculations. Increasing the value of A provides a maximal shift of -0.005 



11 



AE off (meV) 
0.020 ■ 

0.015 

0.010 

0.005 

-0.005 



200 



FIG. 2: (Color online) The energy shift AE a s as a function of the parameter A, using Eq. (39) 



meV, but further increases leads to a positive shift in the energy, due to the dominance of the 
second order terms 5Z 1)2 for large values of A. A positive shift in energy is of the wrong sign 
to explain the proton radius puzzle. 

The requirements of Birse & McGovern [28] can be satisfied by using Eq. (40) with n = 3. 
The use of such a function in calculations of the Lamb shift requires even larger values of A 
to explain the proton radius puzzle. The use of our limit A = 7 leads again to a very small 
increase of the Lamb shift: 0.001 meV. 



V. DISCUSSION 

The principal result we have is that quasi-elastic electron scattering places significant lim- 
its on the off-shell dependence of the nucleon electromagnetic vertex function. While it is 
possible to construct gauge-invariant models of the off-shell behavior that are consistent with 
known features of the virtual photon-proton Compton scattering amplitude, these models are 
incapable of resolving the proton radius puzzle without causing dramatic effects in nuclear 
quasi-elastic scattering in disagreement with observed data. 
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